Analysis of Bifurcations in a Power System Model with Excitation 

Limits 



Rajesh G. Kavasseri and K. R. Padiyar 
Deparment of Electrical Engineering 
Indian Institute of Science, Bangalore, India. 
Email: rajesh.kavasseri@ndsu.nodak.edu 



Abstract 

This paper studies bifurcations in a three node power system when excitation limits are considered. This 
is done by approximating the limiter by a smooth function to facilitate bifurcation analysis. Spectacular 
qualitative changes in the system behavior induced by the limiter are illustrated by two case studies. Period 
doubling bifurcations and multiple attractors are shown to result due to the limiter. Detailed numerical 
simulations are presented to verify the results and illustrate the nature of the attractors and solutions 
involved. 

1 Introduction 

Chaos in simple power system models has been studied extensively in recent papers. In Abed et al., [1993], Tan 
et al., [1993], bifurcations and chaos in a three node power system with a dynamic load model were studied 
using a classical model for the generator. In Rajesh & Padiyar [1999], the authors studied dynamic bifurcations 
in a similar system and reported the existence of chaos even with detailed models. However, in Rajesh & 
Padiyar [1999], it was observed that the field voltage assumed unrealistic values at the onset of chaos owing to 
the unmodeled effect of excitation limits. Though a limiter is fairly easy to model for simulation purposes, the 
effect of a limiter on dynamic bifurcations has been poorly understood because bifurcation analysis demands 
smoothness of the functions describing the model. Limit induced chaotic behavior in a Single Machine Infinite 
Bus system was studied in Ji and Venkatasubramanian [1996] by extensive numerical simulations. In this paper, 
we approximate the limiter by a smooth function to facilitate bifurcation analysis and study the changes which 
arise on it's consideration. The rest of the paper is organized as follows. Section 2 deals with the modeling of 
the system along with the limiter. Section 3 presents the results of a bifurcation analysis. Section 4 contains 
the discussions and Sec. 5, the conclusions. 

2 System Modeling 

The system as considered in Rajesh & Padiyar [1999] is shown in Fig.[T] By a suitable choice of line impedances, 
we might regard the system as one of a generator supplying power to a local load which in turn is connected 
to a remote system modeled as an infinite bus. For the general reader's convenience, a brief explanation of 
the terms d-q and D-Q axis is provided here. The modeling and analysis of three phase synchronous machines 
is complicated by the fact that the basic machine equations are time varying. This is circumvented by the 
use of Park's transformation which transforms the time varying machine equations in to a time invariant set. 
The three phase stator quantities (like voltage, current and flux), when transformed in to Park's frame yield 
the corresponding d-q-o variables. When a generator is described in the d-q frame, then naturally the external 
network connected to it should also be described in the same reference frame. However, the non-uniqueness 
of Park's transformation (each generator has it's own d-q components) prevents us from doing so. In order 
to transform the entire network using a single tranformation with reference to a common reference frame, the 
Kron's transformation where the variables are denoted by D-Q-0 are used. For a complete , detailed and clear 
exposition of these concepts in power system modeling, the reader is refered to Padiyar [1996]. 
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Figure 1: The System 



2.1 Generator Model 
Rotor Equations 

The rotor mechanical equations for the generator as given by the swing equations are, 

5 = LO B S m 

~\- P m Pg 



2H 
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where d is the damping factor in per unit, lob is the system frequency in rad/s, P m is the input power to the 
generator and s m , the generator slip defined by 



LJ — LOB 
L0 B 
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Two electrical circuits are considered on the rotor, the field winding on the d-axis and one damper winding on 
the q-axis. The resulting equations are, 
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The power delivered by the generator P g can be expressed as 

P g = E' q i q + E' d i d + (x'd - x' q )i d i q 

Stator Equations 

Neglecting stator transients and the stator resistance, we have the following algebraic equations 

E' q + x' d i d = v q 
E 'd - x ' q H = v d 



(6) 
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2.2 Excitation System 

The excitation system for the generator is represented by a single time constant high gain AVR and the limiter 
as shown in Fig. 

The equation for this excitation system is given by 

• = -E fdx + K A (V ref -V t ) 
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Figure 2: Excitation System 



E fd = E fdx if Efi n < E fdx < Efr (10) 

rpmin • £ tp ^ rpmin 

~ &fd l J &fdx < &f d 

rpmax ■ r tp - , rpmax 

= E fd if E fdx > E fd 

The limiter shown in Fig. [5] and defined by Eq. 10 is a soft or windup limiter. This limiter model cannot be 
directly used for bifurcation studies. An approximate model where the limiter is described by a smooth function 
is given below (see Fig. |3J. Here, we consider symmetric limits i.e. I-EJ^I = \EJ d n \ — Ef d i 

*2E 

Ef d = fum(E f dx) = — —tarT 1 {aEf dx exp(bE 2 fdx )) (11) 



Figure 3: The function representing the limiter 

Remarks 

Such an approximation amounts to perturbing the vector field slightly and hence the equilibrium structure of 
the system will also be slightly perturbed. So in our studies, the focus will be on how the limiter influences 
non-stationary solutions and their bifurcations. 

2.3 Load Model 

A dynamic load model as in Abed [1993] is used along with a constant power load (Pi d , Qi d ) in parallel with it. 
Thus, the real and reactive load powers are specified by the following equations. 

P = Pid + Pa + piS L + P2 V L + vzVl (12) 
Q = Qid + Qo + qih + q2V L +q 3 V* (13) 

2.4 Network Model 

With the notation defined in Fig. 1, we can write the network equation in the D-Q reference frame as, 

E b + ^r=V t (14) 
Y3 
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Further. 



i = (iq+ jid)e j 



V t = (v q + jv d )e jS 
= h + h 
Y = Y14> = Yi + Y 3 



From Eqs., l(T4*|) to (fTHf> we can write, 



(v q + jv d ) 



Y 



where A x = E b Y^e-^ 5+ ^-^ + Y 1 V L e^ SL - s -^+'f' 1 '> 

and Bx = (iq+jidje-S* 
Defining, 



a = E b Y 3 cos(6 + 4> - fa) + Y x V L cos{5 L - 8 ~4> + fa) 
b = -E b Y 3 sin(6 + (f) ~ fa) + YtVLsin{5t -8-<j) + fa) 



(15) 



(16) 
(17) 
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permits us to write, 

i q cos{4>) + idsin(4>) = Yv q — a 
idCos{4>) — i q sin{4>) = Yvd — b 

2.5 Derivation of the System Model 

Substituting for Vd and v q from the stator algebraic equations (7) and (8), we have, 



(22) 
(23) 
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where Y a = (YE' q - a) 
and Y b = (YE' d - b) 

From Eq. (24), we can solve for the currents id, i q and subsequently solve for Vd and v q from the stator 
algebraic equations. Further, from Eq. l|16f) we get, 

6 = 8 + tan- 1 { — ) 

V a 



(25) 
(26) 



Defining, 



n = s L - 6 ~ fa 

?*2 = 8l - fa 



(27) 
(28) 



the power balance equation at bus 2 can be written as, 

P = VtVLYxcosin) - VlY x cos{fa) + E b V L Y 2 cos(r 2 ) 

-VlY 2 cos(fa) 

Q = VtVLYxsinin) + VlY lS in{fa) + E b V L Y 2 sin(r 2 ) 

+VlY 2 sin(fa) 



(29) 
(30) 
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Substituting from Eqs., (I59T3TI in Eqs., (1-5) and (tTSITTfl) . we get 



x = f(x,A) (31) 

where x = [ 5 s m E' q E' d Efdx 3l Vl ] T and A is a bifurcation parameter. As a simplification, we 
shall also consider the system described the One Axis Model for the generator as the effect of the limiter on this 
case is interesting in itself. For this, we neglect the damper winding on the q-axis and in terms of modeling, 
this is done by omitting E' d as a state variable and substituting 

E' d = -(x q -x' q )i q (32) 

in Eqs., (6) and (8). The state space structure remains the same, with the dimension being one less that the 
previous system. In this case, we have x = [ S s m E' Efdx &l Vl ] • 



3 Bifurcations 

In this section, we illustrate the qualitative differences which arise on consideration of the limiter by studying 

bifurcations in the associated systems with AUT097 (Doedel [1997]) a continuation and bifurcation software 

for ordinary differential equations. The generator input power (P m ) is a very important parameter in practical 

power systems operation. This is the parameter which is adjusted or varied by the power system operators 

(utility) to track the changes and variations in the system load (power demand) so as to maintain a stable 

operating condition. We hence, consider P m i.e. the input power to the generator as the bifurcation parameter. 

To describe the types of bifurcations, we shall use the following notations. 

SNB: Saddle Node Bifurcation 

HB: Hopf Bifurcation 

CFB: Cyclic Fold Bifurcation 

TR : Torus Bifurcation 

PDB : Period Doubling Bifurcation 

In all the bifurcation diagrams the state variable Efdx is plotted against the bifurcation parameter. In the case 
of periodic solutions, we use the maximum value of the variable which is indicated by the circles. Filled circles 
refer to stable solutions and the unfilled ones, to unstable solutions. 

3.1 One Axis Model 
Without limiter 

From Fig. 01 we note that the stationary solutions undergo four bifurcations labeled as HB 1 , HB 2 , HB 3 and 
SNB 4 . For A < A 1 , the equilibrium point is stable, but as A is increased, the stationary point loses its stability 
at A= A 1 through HB 1 . With a further increase in A, the stationary point gains stability through HB 2 , i.e. A= 
A 2 . It remains stable until A= A 3 , where stability is lost through HB 3 . Further, SNB 4 does not influence the 
stability of the stationary point. Next, we focus on the family of periodic solutions emerging from HB 1 . Since 
HB 1 is supercritical, it gives birth to a family of stable periodic solutions indicated by the filled circles. This 
periodic solution loses its stability at TR 5 and with a further increase in A, gains it back through TR 6 and 
remains stable until TR 7 . Further on, there is no qualitative change in its behavior with TR 7 , CFB 8 and TR 9 . 
Next, we find that the branch emerging on continuation of HB 2 is the same as that from HB 1 . On continuation 
of HB 3 , we find a family of unstable periodic solutions which gain stability through CFB 10 . This stable periodic 
solution encounters a PDB 11 on continuation of which, we find a period doubling cascade accumulating at a 
critical value of A c =0.931. (which is not shown here) which definitely suggests the onset of chaos. However, 
what is of interest here, is the behavior of the system after TR 5 . It is clear that a torus bifurcation results in 
the emergence of quasi-periodic solutions. This is verified by simulation as shown in Fig. El which shows the 
quasi-periodic attractor. The bifurcation points are summarized in Table ^ 
With Limiter 

From Fig. El we observe that the stable operating point loses its stability with HB 1 , regains it at HB 2 and loses 
it back at HB 3 before encountering SNB 4 which is similar to the case without limiter (see Fig.^J. Note that in 
Fig- El for the static bifurcations HB 1 - SNB 4 , Efdx < EJ 1 ^ and hence we expect that these bifurcations should 
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Table 1: Bifurcation Points (see Fig. 
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Figure 4: A = P m , One Axis Model 
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Figure 5: The Quasi-periodic trajectory when P m — 0.83 



occur at the same values even with the limiter. However, this is not the case as seen from Table [21 because of 
the approximation which shifts the equilibrium structure as mentioned before. HB 2 and HB 3 occur very closely 
and hence cannot be distinguished in Fig. On continuation of HB 1 which is supercritical, we find that the 
stable periodic solutions do not undergo any bifurcation. HB 2 is also supercritical and its continuation yields 
the same stable periodic set obtained on continuation of HB 1 . HB 3 is sub-critical and its continuation which 
yields CFB 5 where stability is gained for a while before CFB 6 is however, not shown here. Fig. shows the 
time domain plot of the load bus voltage for A = 0.86. The bifurcations are summarized in Table El 



Table 2: Bifurcation Points (see Fig. 
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Figure 6: A = P m , One Axis Model with limiter, continuation of HB 1 and HB 2 
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Figure 7: Sustained oscillations of load bus voltage with time when A = 0.86 with the approximate limiter 



3.2 Two Axis model 
Without limiter 

We let A = P m with reference to Fig. [SJ The stationary point undergoes two bifurcations, HB 1 where it loses 
its stability and SNB 2 which does not influence the stability further. HB 1 is a supercritical bifurcation and 
the family of stable periodic solutions from it undergo a period doubling cascade starting with the PDB 1 , 
accumulating at a critical value of A c = 1.315. The chaotic attractor at A c is shown in Fig. [§] which confirms 
the chaotic nature. The bifurcation points are summarized in Table 
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Table 3: Bifurcation Points (see Fig. [HJ 



Point 
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PDB 3 


PDB 4 
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1.2281 


1.9607 


1.311 


1.314 
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Figure 8: A = P m , Two Axis Model 




Figure 9: The Chaotic trajectory when P m = 1.315 



With Limiter 

From Fig. El we observe that the stable operating point loses stability at HB 1 and then encounters SNB 2 
(which is not shown here). On continuation of HB 1 , which is sub-critical, we find that the unstable periodic 
solution stabilizes with CFB 3 . This stable periodic solution undergoes a period doubling cascade initiated at 
PDB 4 . In Fig. El we also show the period doubled solution and its subsequent bifurcation PDB 5 . By numerical 
simulations, considering both the exact and the function approximation of the limiter, we verify that at A = 1.3, 
the system behavior is chaotic. The time domain plots are shown in Figs ^2 an d El The chaotic attractor 
subject to limits is shown in Fig 1131 The bifurcations are summarized in Table ^ 



Table 4: Bifurcation Points (see Fig. HT]|) 



Point 


HB 1 


SNB 2 


CFB 3 


PDB 4 


PDB 5 


A 


1.2729 


1.923 


1.2557 


1.282 


1.2912 




Figure 10: A = P m , Two Axis Model with limiter 




4 Discussions 



The case studies with two different models were considered solely for illustrating the effect of the limiter on 
bifurcations in the system which is interesting in it's own right. The qualitative difference in system dynamics 
owing to modeling is however not discussed here (see Rajesh and Padiyar [1999] for a discussion). Another 
aspect worth mentioning is the differences in the bifurcation diagrams in this paper from those in the references. 
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Abed et al consider a simplified generator model (classical model) in which the excitation system is entirely 
absent and use a slightly different system for the bifurcation studies. In Ji and Venkatasubramanian [1996], a 
Single Machine Infinite Bus (SMIB) system (which is different from that considered in this paper) wherein the 
load model is absent, is studied. This paper however, focusses mainly on studying bifurcations and changes 
which arise on the consideration of excitation limits. When the One axis model is considered without the 
limiter, we observe stable quasi-periodic trajectories resulting from a TR bifurcation. However, with the limiter, 
we do not observe any bifurcations of periodic solutions with the result that the entire branch from HB 1 to HB 2 
in Fig. is stable. When the Two axis model is considered without the limiter, we observe chaotic trajectories 
due to PDBs, which, with the limiter still occur. However, we observe in this case that the system has multiple 
attractors (see Fig. 110(1 namely, a stable equilibrium point and a stable periodic solution. Further, we observe 
that the PDBs in this case occur very close to the boundary of stable fixed point operation. This means that if 
the system operates close to boundary of stable fixed point operation, and suffers a disturbance with the post- 
disturbance initial condition belonging to the chaotic region, the system can be easily pushed to the chaotic 
region. Another interesting aspect seen by comparing Fig. [H] and Fig. 1101 is that stable equilibrium points close 
to the boundary of stable fixed point operation are surrounded by unstable limit cycles which suggests that the 
region of attraction for the equilibrium points shrinks in the presence of limits. 



5 Conclusions 

An attempt has been made to analyze bifurcations in the presence of a limiter by approximating the limiter 
by a smooth function. It is seen that this methodology provides good insight in to studying bifurcations in a 
system with a soft limiter. The observations in the case studies illustrate in general that, the limiter is capable 
of inducing spectacular qualitative changes in the system. Developing a formal theory for bifurcations and 
analyzing the global system dynamics in the presence of limits in the system would be a challenging task for 
further research in this area. 
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Appendix A 

• Network parameters. 

Fi = 4.9752, Y 2 = 1.6584, Y 3 = 0, fa = fa = fa = -1.4711, E b = 1.0 

• Generator parameters. 

x d = 1.79, x q = 1.71, T' do = 4.3, T' qo = 0.85, x' d = 0.169, 
x' q = 0.23, H = 2.894, w b = 377, d = 0.05, E m = 1.0 

• Load parameters. 

P o = 0.4, ,Q = 0.8 ,pi=0.24 = -0.02 ,p 2 = 1.7 ,q 2 = -1.866 
p 3 = 0.2 ,© = 1.4 

• AVR constants 

K A = 200 , T A = 0.05 

• Limiter constants a = 0.23 , b = 0.1058 , Ef d ax = 5 , Ef d m = -5 
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